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Abstract. We consider the Anderson model on the multi-dimensional cubic 
lattice and prove a positive lower bound on the density of states under cer- 
tain conditions. For example, if the random variables are independently and 
identically distributed and the probability measure has a bounded Lcbcsguc 
density with compact support, and if this density is essentially bounded away 
from zero on its support, then we prove that the density of states is strictly 
positive for Lebesgue-almost every energy in the deterministic spectrum. 



Wegner's estimate, originally formulated in [W] for the Anderson model on the 
lattice Z d , is one of the celebrated tools in the theory of random Schrodinger op- 
erators, see e.g. the recent reviews [V] IKMej or [CHK for the latest developments 
in the case of continuum random Schrodinger operators. In its strongest form for 
lattice models, a Wegner estimate provides Lipschitz continuity of the integrated 
density of states N(E). In particular, this implies that the Lebesgue derivative of 
N(E), the density of states n(E), exists as a function which is essentially bounded 
from above. In addition to the upper bound for the density of states, Wegner also 
presented an argument for a strictly positive lower bound for the density of states 
of the Anderson model in his original article [W]. Although insightful, Wegner's 
argument is not complete as his nonzero lower bound vanishes in the macroscopic 
limit. 

In this note, we give a mathematical proof of a positive lower bound for the 
density of states of the Anderson model. For many years, efforts have been concen- 
trated on Wegner's upper bound because this is essential for the continuity of the 
integrated density of states, the existence of the density of states, and for Anderson 
localization. Consequently, Wegner's idea to obtain a lower bound seems to have 
remained fairly unnoticed. The lower bound is, however, essential for Minami's 
proof that the energy level statistics for energies in the strong localization regime is 
Poissonian [Mm]. Minami fixes an energy E in the region of complete localization 
at which the fractional moment bounds of Aizenman and Molchanov [AMj hold. 
He assumes that n(E) > 0. Minami then proves that the rescaled local eigenvalue 
level spacing measure &Hl{x) = J2j &(L d {ej{L) — E) — x) dx, for the finite-volume 
Hamiltonian (see below for the definition) with eigenvalues Sj (L) , converges in the 
macroscopic limit to a Poisson distribution with density given by n(E). Here, we 
prove the positivity of the density of states at almost every energy in the deter- 
ministic spectrum. We mention that Molchanov [Molj studied the same question of 
energy-level statistics for the one-dimensional Russian school model. In his paper, 
he also proves the positivity of the density of states for that model using completely 
different methods. 
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The Anderson model is given by the discrete random Schrodinger operator H 
on a probability space (0,P) whose realizations := C + lu £ f], act as 

(1) {HW<p)(x) = (£<p)(x)+wMx) 

for all x £ 7L d on a dense domain of ip £ £ 2 (Z d ). Here, the discrete Laplacian £ 
is defined as (Cip)(x) :— J2 y ez d -\x-y\=i fiv)^ an< ^ nas purely absolutely continu- 
ous spectrum cr(£) = [—2d, 2d]. The random potential V^°> consists of a family 
{ U} x}xez d 01 independent, identically distributed real-valued random variables on 

n. 

The Schrodinger operator H of the Anderson model is known [CL1 IPF] to 
be almost surely essentially self- adjoint on the dense subspace {tp £ £ 2 (Z d ) : 
supp tp compact}. Moreover, H is ergodic with respect to lattice translations. To 
define the integrated density of states, we consider finite volumes A C Z d and the 
Dirichlet restriction Hjf' := Ca + V\ of H to the finite-dimensional Hilbert space 
£ 2 (A), where is the restriction of to A and 

(2) (C AV )(x):= ]T tp{y)+tp{x)( J2 H 

y£A:\x~ y\ = l \y(A: |x— y | = 1 / 

for all x £ A and all (p £ ^ 2 (A). Note that the rightmost term in the above 
definition of the Dirichlet Laplacian C\ ensures the Dirichlet-decoupling estimate 
£aiuA 2 ^ £-Ai © £a 2 > see also [51 IKiMii] . We write trA for the trace on £ 2 (A) and 
let Xb stand for the indicator function of some set Bel. Then ergodicity implies 
that the integrated density of states £g1h N(E) is given by the non-random 
limit 

^tr A (X]-^^)) 

along a sequence of expanding cubes A C Z d . Equation ([3|) holds for all E £ R that 
are continuity points of N(E) and all co £ O except for a P-null set, which can be 
chosen uniformly with respect to the aforementioned values of E. 

If the single-site distribution of, say, loq happens to be absolutely continuous with 
respect to Lebesgue measure and if its Lebesgue density p satisfies the additional 
assumption 

(4) p max := ess sup {p(w)} < oo, 

wes. 

then, according to Wegner [W], the integrated density of states N(E) is Lipschitz 
continuous, hence absolutely continuous and the Lebesgue derivative of N(E), the 
density of states, 

(5) £ £ 1 — > n(E) := dN(E) /dE, 

obeys the estimate n(E) ^ p m ax for Lebesgue-almost all E £ M. Another con- 
sequence of ergodicty of the Schrodinger operator H is that there is a closed set 
Scl such that a(H^) — £ with probability one. This set, called the determin- 
istic spectrum of H, is given by £ = [—2d, 2d] + supp p, for the model in JTJ). 
We will not assume Q for the validity of the lower bound for n(E). 

Theorem 1. Let H be the random Schrodinger operator {T]) of the Anderson model. 
Assume that the single-site distribution of u>q is absolutely continuous with respect 



(3) N(E) = lim 
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to Lebesgue measure and that its Lebesgue density p is essentially bounded away 
from zero on some interval [W_, W+] in the sense that 

(6) p mi n := essinf {p(w)} > 

w£[W-,W+\ 

for some —oo < W— < W+ < oo. Assume further that the integrated density of 
states N(E) is an absolutely continuous function with Lebesgue derivative n(E) as 
in ([5]). Then, for every 8 > (small enough) there exists a strictly positive constant 
Cs such that 

(7) n(E) ^C s >0, 

for Lebesgue- almost all E £ [-2d + W- + 8, 2d + W + - 8} . 

Remark 1. The theorem can be generalized in a straightforward manner to in- 
corporate general bounded, self-adjoint and ^-translation invariant unperturbed 
operators H n instead of C. In this case, ([7]) holds for all E £ o-(H )+[W-+5, W+—S\. 

Remark 2. The lower bound constant Cs > can be expressed in terms of the 
integrated density of states (E) for the unperturbed operator C. For any 8 > 
small enough, we cover the interval [—2d + W— + 8, 2d + W+ — 8] by finitely many 
intervals Ij = [Ej — S, Ej + S] of width 28 and centered at Ej . We can take Cs to 
be 

iV (0) (Ej -W-- 28) - iV (0) (Ej -W+ + 28) j , 

where the positive constants Oel > are defined in the proof of Theorem [TJ The 
difference of the integrated densities of states for C on the right of (JSJ) is strictly 
positive, see the end of the proof of Theorem [1] 

In Theorem [TJ we do not require that p is essentially bounded as in ([5]), that the 
support is bounded, nor that [WL,W+] is the entire support of p. However, if wc 
add the latter two hypotheses, we obtain the following special case of Theorem [TJ 

Corollary 1. Under the hypotheses of Theorem^ and the additional condition that 
^ holds on the entire support of p, i.e. p(w) — for almost every w £ R\[W_, W+], 
then the density of states is strictly positive Lebesgue- almost everywhere on S. 

Remark 3. After we completed this work, we learned that a similar result is con- 
tained in Frank Jeske's unpublished PhD-thesis [JJ, which was supervised by Werner 
Kirsch. We thank Ivan Veselic for informing us about the existence of [JJ. 

Remark 4. The question arises naturally whether a similar lower bound for the 
density of states n does also hold in the case of continuum random Schrodinger 
operators, that is, Schrodinger operators on L 2 (R d ). For d = 1 and for alloy-type 
random potentials with suitably well-behaved single-site potentials, the answer is 
affirmative. The argument proceeds as in the discrete case with some obvious mod- 
ifications that are well known from proofs of upper Wegner estimates for continuum 
models. The key point is that the finite-rank-perturbation argument, which allows 
us to proceed from (|15[) to (|16[) below is still valid in the one-dimensional continuum 
case. Indeed, suppose we have two Schrodinger operators on an interval that differ 
only by a boundary condition (Dirichlet versus none, say) that is imposed at an 
interior point of the interval. Then it is well known from, e.g., the theory of point 
interactions that these two Schrodinger operators differ by a rank-2 perturbation. 



(8) Cs = mm — 
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However, for d ^ 2, different boundary conditions along a finite hypersurface S give 
rise to an infinite-rank perturbation. Thus, in the case d ^ 2 one needs an alterna- 
tive argument why different boundary conditions along S for Schrodinger operators 
in a finite volume A (with S in the interior of A) lead to eigenvalue counting func- 
tions that differ by a term proportional to the area of S. Furthermore, this error 
term would be required to remain bounded as A f M. d . But this is a delicate issue 
in view of [KT1IK21. 



Proof of Theorem[I[ The proof follows Wegner's arguments [W], except that we 
introduce a partition of the finite volume into cubes of large but fixed size. This 
allows us to get a nontrivial result in the macroscopic limit L — > oo, a problem 
seemingly ignored in [Wj . A similar partitioning strategy was used in the proof of a 
(upper) Wegner estimate for continuum random Schrodinger operators by spectral 
averaging [CH] , see also |FHLM| IHLMWj for the case of Gaussian or other types 
of unbounded random potentials. 

1. Let Ei, E 2 & M. such that E2 — E\ > e for some e > 0. We consider a sequence 
of expanding cubes A^ in Z d with volume \Al\ = L d . Finally, we pick a smooth, 
monotone increasing switch function f E € C 1 (R) such that / £ (A) = for all A ^ 
and ,f e (A) = 1 for all A ^ e. We let E denote the expectation associated with the 
probability measure . 



N(E 2 ) - N(Ei) ^ lim 

L — >oo 



and we write F e ,L (A, lu) 
1 



— trA £ / S (A — Hj^). Then we have 
-H AL )-f e (E 1+ e-H 



(9) 




(fe(E 2 



Ai 



)) 



——; I dAE 



The quantity -F ei L(A, w) depends on A and w only through the differences {u> x — 
\} xe A L , and it is a monotone decreasing function in each of those differences. We 
partition the cube A^ into (L/£) d smaller cubes Fj of the same (fixed) volume I . 
We consider only those big cubes A^ for which such a partition is possible. We will 
take L — > 00, and I large but finite. Therefore we get 



(10) 



d_ 

ox 



F s<L (\,u) 



(L/£) d 

E E 

3=1 xeTj 



d 
duj x 



for all A G R and all lu G fl 

2. We conclude from (O and Gil that 



(11) N(E 2 )-N(E 1 )^pi 



lim 



{L/C) 



(L/l) d 

, E *n 

3=1 



E 2 



dA Gj(X, 



with 

(12) G,(\,u„) 



Here T c - := lj d \ Tj denotes the complement of Tj, and (in slight abuse of notation) 
we have written lo =: (wr^,, ^r^), where LOVj '■= {io x ) xe r j - The partial disorder 
average Er c in (fT2|) extends only over the coupling constants wr c • 
3. Following Wegner [W], we are going to perform a change of variables in 
(|12| from LOTj to 77: we fix an arbitrary point Xj G Tj and set r] Xj :— to Xj and 
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rj y := ujy — u> Xj for all y G Tj \ {xj}. The Jacobian associated with this change of 
variables is 1, whence 

G M^) = y d f d Vxj I ( n d %) 

I J[w.,w + ] Jiw--v* r w + -v* j ]* d --\ e r,-\ {xj } ' 

(13) x (-^^(A^faWj)). 

Now, fix 5 G]0, (Wlf — M / _)/4[. One obtains a lower bound for (fT3"|) by restricting 
first the integration over r) x . to [W_ + 5/2, W+ — (5/2] and then restricting the 
integration over r\ v to [—6/2, 6/2], for all y G Tj \ {xj}. This gives 

« ■/[-«/2,«/2|^- 1 \ 6 r7\{* J } 71 

(14) -F E , L (A,(w r ,(r ? +),cJr|)) 

with ?7 ± := (W± T <V 2 > (%) !y er J \{^})- Note that in 03} one has (ur^r)')) x < 
VF- + J and (tUTj (?7 + )) x 5^ W + — 6 for all x G Ij . Since F £j l is a decreasing function 
in each uj x , we arrive at 

(15) G^urc) > ^^(A.^.wr.)) - F e , L (A, (w+ >Wr .)) 



with spatially constant couplings o; r . := (W± =F <5)ser : , inside the small cube Tj. 
4. Next, we will use a Dirichlet decoupling of the small cube Tj. In the first 
(i.e. the positive) term on the right-hand side of (fT5"|) . this can be done straight 
away, because H\ L ^ £Zrv © H\ L \^.. To do the replacement in the second (i.e. 
the negative) term, one has to take into account the error that arises from intro- 
ducing the additional Dirichlet boundary condition along dTj \ dA^. But this is 
a perturbation of rank 0{l d ~ l ) and it is independent of the coupling constants. 
Furthermore, recall that ^ f e ^ 1. Thus there is a constant D G]0,oo[, which 
depends only on d, such that 

x£ d — l 

G j (A, wr«) ^ — {tr r , [f £ (X -W- - 6 - C Tj ) - f c (X-W+ + 6 - C T] )] 
-Di d ~ 1 \ 

xt d — i 

> — { tr r, [fe(Ei -W--5-C Fi )- f s (E 2 -W++6-C T] )] 

(16) -Dl d -^ 

for all A G \Ei,E2\ and all Wr?- The contributions from H^^^ have canceled, so 
the right side in (fT6| is independent of L. Inserting (fl6|) into (jlljl and taking the 
limit e | 0, we arrive at the estimate 

(17) m^mi , (i£ft (Kt{EuE2) _ D/t) 

for the difference quotient of the integrated density of states of H. The lower bound 
in (|17|) is expressed in terms of the difference 

(18) K t {E x , F 2 ) := N$ {E x - W_ - <J) - JVjjJ {E 2 -W++6) 
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of the Dirichlet finite-volume approximation (A) := £~ d ti'A f X]- 00t x\(£A e ) for 
the integrated density of states of the free Laplacian C. 

5. By hypothesis we know that N(E) is absolutely continuous with respect to 
Lebesgue measure. Hence, we can take the monotone limit E2 | E\ =: E in (TT7|) 
and obtain 

(19) n(E) > {dPm ^ Y (Kt(E,E) - D/l) 

for Lebesgue-almost all £eR. We now fix E S K, and observe that Kg(E, E) ^ 
Ki(E -5,E + 5)=: K e (E Q ), for all E E [E Q - 5, B + 5]. We next note that 

(20) K{E ) := lim i^(£J ) = iV (0) {E -W-- 25) - N {0) (E -W+ + 25) 

£—>oc 

exists, where iV(°)(A) := lim^oo (A) = (5q, X]_cxj,A] (£)8o). It is important to 
observe that E - W- - 28 > E -W + + 25, since < 8 < (W + - W_)/4, and 
that for all E e] - 2d + W_ + 25, 2d + W + - 2S[, we have -2d < E - - 25 < 
2d+(W+-W-)-A5, and -2d- [(W+ - W-) -A5] < E -W+ + 25 < 2d. Specifically, 
we have — 2d < E — W- — 25 and if E — W- — 25 > 2d, then the other energy 
satisfies Eq — W+ + 25 < 2d. Consequently, K{Eq) is strictly positive for Eq on the 
specified range since 

(21) AT ( °)(A 2 ) -iVW(Ai) > 0, 

whenever Ai < A2 , and at least one of the Xj 's lies in the interior of <j(£) = [—2d, 2d] . 
Thus, there exists a finite length £e such that 

(22) K tEo (E, E) - D/l Eo > K(E a )/2 for all E G [Eq — 5, E + 5}. 

The theorem follows from (fT9)) , ([22]) and by covering the interval ] — 2d + W- + 
5, 2d + W+ — 5[ by a finite number of small intervals of length 25. □ 

Acknowledgment. The authors thank Frangois Germinet for the kind hospitality at 
the Universite de Cergy-Pontoise, and Abel Klein for several discussions. 
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